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Abstract 

Unsteady interactions of strong concentrated 
vortices, distributed gusts, and sharp-edged gusts 
with stationary airfoils have been analyzed in 
two-dimensional transonic flow. A simple and 
efficient method for introducing such vortical 
disturbances has been implemented in numerical 
codes that range from inviscid transonic small- 
disturbance to thin-layer Navier Stokes. The 
numerical results demonstrate the large distor- 
tions in 'the overall flow field and in the surface 
air loads that are produced by various vortical 
interactions. The results of the different codes 
are in excellent qualitative agreement, but, as 
might be expected, the transonic small-disturbance 
calculations are deficient in the important region 
near the leading edge. 


Introduction 


Vortical disturbances can have important 
effects on the airloads and the aeroacoustics of 
a variety of aerodynamic devices, such as maneu- 
vering aircraft and missiles, helicopter rotor 
blades, and turbomachinery. Although the numeri- 
cal analysis of strong vortical flows has received, 
less attention than the treatment of shock waves 
over the past decade, compressible flow fields 
with embedded regions of distributed or concen- 
trated vorticity appear certain to receive 
increasing attention in the future. 

Several methods have already begun to appear 
in the literature for treating concentrated vor- 
ticity in finite-difference computations. Direct 
approaches, such as Rizzi and Erickson, 1 Fuiii and 
Kutler, 2 Krause, ^ and Srinivasan and Steger 1 * may 
be classified as vortex-capturing, by analogy with 
the familiar shock-capturing methods that are used 
in many computational fluid dynamics (CFD) codes. 
As in shock capturing, the details of the actual 
phenomenon are spread over several grid points. 
Consequently, the solutions in these high-gradient 
regions are artificially grid-dependent and are 
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susceptible to the effects of numerical dissipa- 
tion. In the case of concentrated-tip vortices, 
for example, numerical dissipation may destroy the 
simulation of the core structure faster than would 
physical dissipation, 11 unless many grid points are 
clustered in the vortical region. 

Two alternate methods for modeling vortex 
flows are illustrated schematically in Fig. 1 . By 
analogy with shock-fitting in transonic flows, 
these methods may be thought of as vortex-fitting. 
The upper part of the figure portrays the method 
introduced by Caradonna et al.,^ and used by 
George and Chang, ^ for example. This technique of 
one or more branch cuts permits concentrated vor- 
tices to be introduced into potential-flow 
formulations . 


The prescribed-vortex method, which is the 
subject of this paper, is shown in the lower half 

of Fig. 1. The essential ideas are due to 
7 P 

Steinhoff' and Buning and Steger, and the tech- 

Q-1 1 

nique has been used by the authors 7 for poten- 
tial, Euler, and thin-layer Navier-Stokes analyses 
of two-dimensional, unsteady airfoil-vortex inter- 
actions. In these applications, the structure of 
the vortex is prescribed, but its path in space 
can be allowed to develop as part of the solution. 
The prescribed-disturbance method is also well 
suited to and has been used for more general vor- 
tical flows, such as gusts with distributed 
vorticity. 1 2 


Many of the interesting practical applica- 
tions, such as blade-vortex interactions of heli- 
copter rotors (e.g., Fig. 2), and gust response of 
aircraft wings, are three-dimensional unsteady 
problems that will require enormous computational 
resources for transonic cases. Therefore, it is 
useful to study the numerical methods, as well as 
the general nature of the phenomena, in two dimen- 
sions, as indicated in the figure. It is also 
important to establish the minimum level of com- 
plexity in the governing equations that will suf- 
fice. In the following sections, we review our 
transonic small-disturbance, Euler, and thin-layer 
Navier-Stokes analyses, especially the two- 
dimensional implementation of the prescribed- 
disturbance approach. Several representative 
examples are given that illustrate the advantages, 
capabilities, and limitations of these different 
formulations, in anticipation of future extensions 
to real, three-dimensional problems. 
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BRANCH CUT 
TO DEFINE 
VORTEX, [0] = r v 


WAKE SHEET 
[0] = r w (y) 



(a) POTENTIAL - FLOW VORTEX FITTING 



r y CALCULATED (FORCE-FREE PATH) 


(b) PRESCRIBED - DISTURBANCE VORTEX FITTING 

Fig. 1 Methods for modeling vortex flows. 


2-D AIRFOIL, A, = 0 



Fig. 2 Two- and three-dimensional unsteady vortex interactions. 


Numerical Formulations tion of the governing equation set. Even though 

7 the governing equation(s) may be nonlinear and 

The essence of the prescribed-disturbance' or indepe ndent solutions are therefore not super- 
perturbation 0 method is simply that each of the p 03 able, we may still decompose the dependent 

dependent flow variables is split into a pre- variables as follows: 

scribed part, which defines (or is consistent 
with) the imposed vortical disturbance, and a 
remaining part, which is obtained from the solu- 
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Le J 

and q v represents the specified vortical distur- 
bance superimposed on a uniform free stream. 

The following subsections summarize the 
development of the governing equations for 
(q " q v )» and present the methods for solving them 
numerically. Although the overall concept is 
essentially the same for the Euler or Navier- 
Stokes formulations as for the transonic small- 
disturbance (potential) formulation, the actual 
derivations and solution techniques are rather 
different; therefore, they are treated separately. 
The complete details of each are given in Ref 3. 10 
and 11 for the former, and in Refs. 9 and 12 for 
the latter. 

Euler and Thin-Layer Navier-Stokes Equations 


♦ u(n x + n y )(u2 + y2) n /2 
L + (y/3)(n x u + n y v)(n x u + n v )J 

Here U and V are the eontravariant velocities 
along the £ and n directions, i.e., 

U - + 5 x u + 5 y v (5a) 

V - n t + n x u + n y v (5b) 

and J is the transformation Jacobian. 1 ^ 

In Eq. (2) q y represents any prescribed 
nonuniformity superimposed on an otherwise uniform 
free stream. In the present investigation, a 
concentrated vortex is considered, either fixed in 
space or moving with the flow. The velocity field 
of the vortex is specified as 


In this part of the investigation, we use the 
unsteady, two-dimensional equations in strong 
conservation-law form, 1 ^ including the Baldwin- 
Lomax 1 ^ algebraic eddy-viscosity model for thin- 
layer turbulent flows. In generalized coordinates 
(5,n,i) and in the perturbation form of Buning and 
Steger,® the equations are 

3 T (q - %) ♦ 3 5 (E - E v ) ♦ 3 n (F - F v > - Re' 1 3^ 


( 2 ) 


where the dependent flow-field vector is 
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and the flux vectors E, F, and S are given by 
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( 3 ) 


2 tt r 2,2 

v(9) = (1 - e r /a ) (6) 

|r| 

where^ T y is the dimensionless vortex strength 
and r(x,t) is the instantaneous radial distance 
from the vortex center. That is, the vortex- 
induced field q y is irrotational outside a vis- 
cous core of radius a. Thus the vector compo- 
nents p, (pu), (pv), and e are determined from 
the inviscid (Euler) components of Eq. (2), as 
discussed in Refs. 10 and 11. 

The usual tangency and no-slip boundary con- 
ditions on the airfoil are applied at each instant 
to the total flow-field vector q. Surface- 
conforming C-type grids are used, generated by an 
algebraic method, ^ with clustering near the lead- 
ing and trailing edges and in the vicinity of the 
upper and lower surface shock waves. Typical 
grids for viscous calculations consist of 1 6 1 to 
181 points in the surface, 5, direction, and 
52 points in the normal, n, direction, extending 
6 to 10 chords in all directions. Forty-five grid 
points in the normal direction are generally used 
for the inviscid cases. 

An implicit, spatially factored algorithm 
with Euler-implicit time-differencing 1 ^ is used to 
solve Eq. (2). The scheme is first-order accurate 
In time and second-order accurate in 3pace. Sec- 
ond-order Implicit smoothing and fourth-order 
explicit dissipation terms are added to improve 
the nonlinear stability. Even so, the nondimen- 
sional time steps generally have to be restricted 
to values on the order of 0.005 during that por- 
tion of the interaction when the vortex is close 
to the airfoil. The CPU solution time on the Ames 
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Cray X-MP computer is approximately 2 sec/time 
step. The inviscid cases considered in this paper 
required approximately 10-15 min each, and the 
viscous calculations about 25-30 min CPU time, 
after the initial, steady solution was obtained 
without the vortex. 


however, the continuity and momentum equations can 
be combined to give 12 


3Q 

3t 
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U > 37 
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(7) 


where 


The effectiveness of the present prescribed- 
disturbance scheme is illustrated in Fig. 3 for 
one of the examples to be discussed in more detail 
in the section on Unsteady Vortex Interactions. 

The figure shows the variation in lift on the 
airfoil as a strong, concentrated vortex passes 
beneath it, with and without the special treatment 
outlined above. The "nonperturbation" scheme 
essentially represents vortex capturing; that is, 
the vortex structure is specified only as an ini- 
tial condition, when the center of the vortex is 
5 chords upstream of the leading edge. Subse- 
quently, this initial disturbance convects through 
the computational domain, and the numerical dissi- 
pation associated with the finite grid spacing 
progressively weakens the gradients and reduces 
the effective vortex strength substantially. This 
numerical error is clearly grid-dependent; how- 
ever, it is completely absent in the prescribed- 
disturbance solutions, which are essentially inde- 
pendent of the grid. 


181 x 45 GRID 


NON-PERTURBATION 

PERTURBATION 



Fig. 3 Effectiveness of the prescribed- 
disturbance method for Euler calculations. 


Transonic Small-Disturbance Potential 
Equation 


Q 


■/(If - s * » « ;) 


• dr - f(t) 


( 8 ) 


dr is an incremental position vector which is 
integrated over a suitable path from some upstream 
reference point, and the symbol q now represents 
the velocity field (u,v). 

We now go a step beyond Eq. (1) and decompose 
this total velocity field into the following com- 
ponents: 1) the uniform free stream, U m ; 2) the 

prescribed disturbance field, q ; and 3) an irro- 
tational perturbation velocity field V<|> produced 
by the airfoil in the presence of the concentrated 
vortex or distributed gust. That is, 

q - 5, + q v + v* (9) 

Then, for example, for a distributed vortical 
(rotational) gust 

Q ’/(sT- ' ^ x 7 x ^v) ’ d? + If (,0a) 


whereas for an irrotational moving vortex 


Q 


♦ 3i 

3t at 


(10b) 


where <f> y * 2irT tan* 1 C(y - yj)/(x - x^)], and 
(Xj,yj) is the instantaneous position of the 
center of the concentrated vortex. 


It should be emphasized again that the decom- 
position of the velocity in Eqs. (1) or (9) does 
not imply linearity. Both the boundary conditions 
and the governing equation for the airfoil distur- 
bance potential, V$, are altered by the introduc- 
tion of q y , and independent solutions are not 
superposable for transonic flows. 


In this case, the approximations are a thin 
airfoil in an inviscid, isentropic fluid, with 
either (1) a concentrated potential (irrotational) 
vortex, or (2) small distributed vortical distur- 
bances, superimposed on an otherwise uniform, 
nearly sonic free stream. As in the above formu- 
lation, these prescribed vortical disturbances are 
assumed to not be altered by the interaction with 
the airfoil. 


The governing equation for <£ is^obtained by 
substituting first the quantity q Q = + q y , and 

then q given by Eq. (9), into Eq. (7). Next, 
the former resulting equation is subtracted from 
the latter. Then the quantity 

/(ar - % x v x q v ) . a? 


The combination of unsteady flow and rota- 
tional disturbances in the free stream requires 
some special attention if potential-flow concepts 
are to be retained. Under the previously men- 
tioned isentropic, small-disturbance assumptions, 


which is common to both, drops out. Even if 
V x q y • 0, the quantity (7$ x V x q ) -dr 
vanishes if dr is aligned with V<t>. Finally, 
the local speed of sound, a, is evaluated using 
the isentropic relations and the usual tran 
sonic small-disturbance scaling laws and 
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approximations, 16 ” 17 coupled with the observation 
that the prescribed disturbances considered in 
this paper produce only second-order variations in 
density, pressure, and temperature. This gives 
the following modified form of the unsteady 
transonic small-disturbance equation (in strong 
conservation form) 


M -»tt + 2M -*xt - h C(s2 + C 2*X )( *X + V 


- 6 u ] + <J> 

v yy 


(ID 


where 


8 2 - 1 - M 2 


c 2 ■ - j Or - Dm! 


Linear theory can be simulated by setting C 2 * 0. 

The boundary condition on the airfoil is 
obtained by simply substituting Eq. (9) into the 
small-disturbance approximation to flow tangency, 
giving 



v v (12) 


As noted above, neither small-amplitude gusts 
nor moving concentrated potential vortices produce 
first-order pressure fluctuations in the free 
stream. Therefore, consistent with our basic 
assumption that the prescribed vortical distur- 
bance remains undistorted as it passes by the 
airfoil, we retain the usual form for the small- 
disturbance pressure coefficient, 

C p = -2(<fr x ♦ * t ) (13) 

Numerical experimentation revealed that the 
most serious deficiencies of the small-disturbance 
formulation are 1 ) a limitation on the strength of 
the vortex and/or its closest proximity to the 
airfoil, and 2) the leading-edge behavior for 
strong interactions with very thin airfoils. The 
latter can be partially overcome with an ad hoc 
leading-edge correction to Eq. (12); viz. 


&). - V x" r 7 28 [(ft + u » h) Y b - v v] (12b) 


as discussed in Ref. 10. This treatment has been 
found to be very satisfactory for thin airfoils, 
such as the NACA 64A006 considered in the section 
on Unsteady Vortex Interactions. However, compar- 
isons with Euler and full-potential solutions show 
that it over-corrects for moderately thick air- 
foils, such as the NACA 0012, and that it is com- 
pletely unsatisfactory for airfoils with blunter 
leading edges. 


Equation (11) is solved using the alternat- 
ing-direction-implicit numerical algorithm of the 
well-exercised transonic code LTRAN2 , 1 7 with all 
of the high-frequency terms included 9 and provi- 
sions added for variable time steps. A rectangu- 


lar cartesian grid (113 or 186 x 97) is used, with 
clustering near the leading and trailing edges. 

The CPU solution times on the Ames Cray X-MP com- 
puter range from approximately 1 to 4 min. 

Comparison of Initial Solutions 

A preview of the similarities and differences 
in representative initial solutions, representing 
the steady-state solutions of airfoils in uniform 
free streams, Fig. 4, is helpful before consider- 
ing the vortical interactions in the section on 
Unsteady Vortex Interactions. In this and subse- 
quent figures, the Reynolds number is 6 x 10 6 for 
the thin-layer Navier-Stokes solution, with a 
turbulent boundary layer assumed for the entire 
airfoil. The leading-edge correction in Eq. (12b) 
was used for the small-disturbance calculations on 
the NACA 64A006 airfoil, but not for the 
NACA 0012. The small-disturbance deficiency in 
the leading-edge region is more evident in the 
results at M # * 0.3 than at higher Mach number. 

For the NACA 0012 airfoil at » 0.8, the 
shock wave is fairly strong, but viscous effects 
are relatively weak. Therefore, the Navier-Stokes 
and Euler solutions are approximately the same, 
but the details near the shock wave are slightly 
different in the small-disturbance solution. On 
the other hand, the viscous effects are stronger 
on the 64A006 airfoil, but the Euler and small- 
disturbance are almost identical. The experimen- 
tal results are from Zwaan. 16 


Example Results — Unsteady Vortex Interactions 

For the sake of brevity, we consider only 
concentrated vortices convecting with the flow 
past symmetrical airfoils at zero angle of attack. 
Therefore, all the differences between the flow 
fields above and below the airfoils are due solely 
to the vortex interaction. Examples of vortices 
fixed in space and of airfoils at incidence are 
given in Refs. 9-1 1 . 

The small-disturbance results in Fig. 5 illu- 
strate many of the essential features of transonic 
airf oil-vortex interactions. For this case of a 
vortex with clockwise circulation passing beneath 
the airfoil, the vortex-induced "downwash" effect 
first appears as an effective negative angle of 
attack. Figs. 5a and 5d. Unlike the linear 
results, however, the transonic effects are much 
larger on the lower surface than they are on the 
upper surface. Also, the actual shape of the 
pressure distributions, the shock-wave strengths 
and positions, and the time-history of the airfoil 
pressures are quite different, even qualitatively, 
from the results for airfoil oscillations, or from 
the results for the sinusoidal gusts considered in 
the section on Airfoil-Gust Interactions. 

As the vortex passes beneath the airfoil, 
Figs. 5b and 5e, the u-component of the vortex- 
induced velocity field also becomes important in 
the nonlinear case. After the vortex passes 
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NACA64A006 AIRFOIL, = 0.85, a = 0° 

- 8 r r r 



Fig. 4 Initial solutions without a vortex. 


behind the trailing edge, Figs. 5c and 5f, its 
induced field becomes and effective "upwash" 
distribution of vertical velocity, and its 
horizontal component becomes relatively unim- 
portant again. Throughout the interaction, non 
linear unsteady lag effects are very important, 


and the return of the flow around the airfoil to 
Its orginal state is extremely slow. The moving- 
vortex interaction is substantially different from 
the quasi-steady response that would be obtained 
by freezing the vortex at different x-positions. 
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Fig. 5 Transonic airfoil-vortex interaction on the NACA 0012 airfoil at 

- 0.8, a - 0, r - 0.20, y - -0.26. 

ou v v« 


Figures 6 and 7 compare the transonic small- 
disturbance, Euler, and thin-layer Navier-Stokes 
results for the NACA 6*JA006 airfoil at M w - 0.85, 
with the same vortex strength and initial position 
as before. This airfoil, which has been used in 
numerous numerical and experimental studies, is 
not only thinner than the NACA 0012 section, but . 
has a significantly smaller leading-edge radius as 
well. Therefore, it might be expected to be more 
sensitive to the vortex-induced downwash than the 
NACA 0012 profile. 

As in the steady results shown in Fig. *1, the 
shock wave in this case is smeared by the pre- 
dicted viscous-inviscid interaction, according to 
the Navier-Stokes calculations. However, no 
boundary-layer separation is induced. Otherwise, 
the predictions of all three methods are similar, 
and the sequence of events is essentially the same 
as those discussed in connection with Fig. 5. 
However, larger fluctuations appear in the pres- 
sure distributions near the leading edge. Without 
the leading-edge correction of Eq. (12b) this 


effect was over-predicted by the small-disturbance 
calculations. 

For the results shown in Figs. 6 and 7, the 
vortex is constrained to convect at U w along a 
straight line, y v - -0.26. Figure 8 shows the 
effect of allowing the vortex to move along a 
force-free path; that is, along a streamline. “The 
Euler calculations show some differences caused by 
this effect as the vortex passes near and beyond 
the trailing edge, but this behavior is not 
evident in the small-disturbance results. 

Instabilities usually develop in the small- 
disturbance calculations when cases with sig- 
nificantly stronger vortex interaction are 
attempted. The Euler and Navier-Stokes calcu- 
lations require a reduction in the nondimensional 
time-step and require special attention to the 
dissipation and smoothing of parameters. However, 
solutions have been obtained 1 ^' 11 for values of 
f y up to twice the values given in Figs. 5-8. 
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1.0 


THIN-LAYER NAVIER-STOKES 


x v = -0.3 x y = 0 I x y = 0.5 I x v = 1.0 



Fig. 6 Airfoil-vortex interaction on the NACA 6HA006 airfoil at 
M - 0.85, a - 0, r - 0.20, y - -0.26 « constant. Re - 6 Mil. 

OO * V "v 


THIN-LAYER 
o NAVIER-STOKES 


EULER 

A SMALL 



Fig. 7 Lift- and pitching-moment variations with 
vortex position on the NACA 64A006 airfoil at 
M a - 0.85, a - 0, r - 0.20, 
y v - -0.26 - constant. 


EULER 


PRESCRIBED PATH 
FORCE-FREE PATH 


ATRAN2 ^ 

o 


PRESCRIBED PATH 
FORCE-FREE PATH 




Fig. 8 Lift- and pitching-moment variations with 
vortex position on the NACA 64A006 airfoil at 
M» - 0.85, a - 0, f - 0.20, y - -0.26. 
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Examples of Airfoil-Gust Interactions 

A number of calculations have been performed 
with the small-disturbance code for transonic 
flows with sinusoidal and sharp-edged gusts, whose 
distributions of vorticlty are small in amplitude, 
but which are not necessarily irrotational . 12 
These correspond to examples that have previously 
appeared in the literature for incompressible or 
subsonic flows. 

Sinusoidal Transverse Cust3 

The first case to be considered is the tran- 
sonic counterpart of the well-known sinusoidal 
traveling-wave gust with a transverse velocity 

v Q - A 1 sin(wt - ojx) (14) 

where x/c - 2x/c - 1 is measured from the mid- 
chord of the airfoil, and where c i3 the chord. 
This case is particularly straightforward within 


the framework of the transonic small-disturbance 
approximation, since the gust velocity does not 
appear explicitly in the governing equation. The 
effect of the gust is felt through the boundary 
condition on the airfoil as an effective time- and 
space-dependent angle of attack, Eq. (12). 

A sketch of this problem and typical results 
for this type of gust are shown in Fig. 9; the 
figure shows computations for both the nonlinear 
and linear mode of Eq. (11) at various "phases” 
during the passage of the gust. This figure may 
be compared with Fig. 5, for which the airfoil, 
Mach number, and mean angle of attack are the 
same. Here the nondimensional frequency param- 
eter, K - ioc/U^, is inversely proportional to the 
wave length X of the gU3t; X - 2-irU^/to - 4irc in 
this particular case. The amplitude of the gust 
corresponds to an effective fluctuating angle of 
attack of ±1°. 


NACA 0012 AIRFOIL, = 0.8, a = 0°, A 1 » 1°, coc/U^ = 0.5 







X 


Fig. 9 Transonic airfoil-gust interaction on the NACA 0012 airfoil at M w - 0.8, a - 0; sinusoidal gust 

with A 1 - 0.017, wc/U^ - 0.5. 
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The main distinction of the linear and non- 
linear solutions in Fig. 9 i3 the presence of a 
shock wave and its unsteady motion as the gust 
passes by the airfoil. The role of moving 3hock 
waves is even more dramatic in the following exam- 
ple, which is a comparison of the conventional 
NACA 0012 airfoil and a modern supercritical 
profile . 


Computed results for the Dornier CAST-7 and 
NACA 0012 airfoils with significant amounts of 
lift are compared in Figs. 10 and 11. Figure 11 
shows the harmonic components of the surface pres- 
sure distributions, referenced to the phase and 
amplitude of the prescribed gust. The nonlinear 
effects are much stronger on the upper surface 
than in the previous example. It is interesting 
to note that upper-surface shock wave is much 
stronger on the NACA 0012 than on the CAST-7 air- 
foil, with correspondingly larger magnitudes of 
pressure fluctuations in the neighborhood of a 
shock wave. However, the shock-wave motion is 
much greater on the CAST-7 airfoil, to the extent 
that the mean C p distribution in Fig. 11 is 
significantly different from the C p distribution 
without the gust, Fig. 10. 


SUPERCRITICAL 
CAST-7 AIRFOIL 
a = 0° 

C, 2: 0.70 
L o 



Fig. 10 Steady Transonic 


and conventional airfoils 


CONVENTIONAL 
NACA 0012 

a = 2.5° 

C, a* 0.70 
u o 



flow on supercritical 
at M m = 0.76. 


The lower surfaces of both airfoils remain 
subcritical throughout most of the cycle of the 
gust interaction. Consequently, both curves of 
mean C p in Fig. 11 are approximately equivalent 
to the respective lower-surface steady values in 
Fig. 10 and there is no higher-harmonic content. 
There are only mild distortions of the first- 
harmonic fluctuating pressures in the vicinity of 
the shock wave, and the overall features approxi- 
mately resemble the linear solutions (not shown). 


CAST-7, <* n - 0, C, - 0.7 

o U Q 

NACA 0012, cx Q = 2.5°, C L - 0.7 

° *-0 




(a) Upper surface. 


CAST-7, a =0, C, = 0.7 
° L o 

NACA 0012, q q = 2.5°, C Lq s 0.7 



(b) Lower surface. 


Sharp-Edged Gusts — Indlcial Response 

Indicial responses to impulsive plunging 
motion, step changes in angle of attack, and 
sharp-edged gusts are often used in linear analy- 
ses to determine the aerodynamic loads for arbi- 
trary body motions or gust fields. Ballhaus and 
Goorjian 1 ^ have discussed how this approach can 
be extended to a limited class of transonic 


Fig. 11 Harmonic components of the pressure dis- 
tributions on supercritical and conventional air- 
foils at M - 0.76; sinusoidal gust with 
A 1 - 0.017, wc/U^ * 0.5. 

oscillating airfoil problems, if the shock-wave 
motion is negligible and the unsteady fluctuations 
are essentially linear perturbations about a non- 
linear steady or mean flow. In this section we 
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discuss the indicial response to a sharp-edged 
gust, Fig. 12, in subsonic and transonic flow. 



V G 


SHARP-EDGED 

GUST 




Fig. 12 Schematic of indicial response for 
plunging motion (above) and for a sharp-edged gust 
(below) . 


Figure 12 illustrates the similarities and 
differences in the response of a thin airfoil to 
two types of indicial inputs; namely, 


plunge: a 


gust: Vq 


{ 0 for t < 0 

A 1 for t > 0 

f 0 for U t < X 

j 

l A, for UJ; > X 


05) 

( 16 ) 


20 

According to linear theory, the sharp-edged gust 
response rises monotonically from zero at t - 0 
to the value C. - 2 ttA_/( 1 - M 2 )' as t -► ®, 
whereas the plunge response starts at itA 1 

for M - 0 and at f)A,/M for M >0. In 
either case, the approach to the steady-state 
value is progressively slower as M w increases. 


Figure 13 shows the nonlinear gust response 
of the NACA 64A006 airfoil as a function of Mach 
number, with a logarithmic abscissa. The critical 


Mach number for this airfoil is approximately O .83 
for o « 0 , and the flow is only weakly transonic 
at M w - 0.80 and a - 1°. Therefore, the lift 
response shown in Fig. 13 is essentially linear 
for M w < 0 . 80 . 

However, the behavior becomes distinctly 
nonlinear at higher Mach numbers. The overshoot 
in lift at M w - 0.9 is due to a complex, non- 
linear motion of the shock wave on the lower sur- 
face, as shown in Fig. 14. The upper-surface 
shock wave moves monotonically rearward from its 
initial position of x » 0.85 and stabilizes at 
the trailing edge at U tt t/c * 25. However, the 
shock wave on the lower surface moves forward 
initially, reaches a minimum value of X s at 
U^t/c » 30, and then undergoes a small oscillation 
in position between U^t/c » 40 and 60 , before 
reaching its final rearward location. This 
nonlinear behavior indicates that it would be a 
mistake to use the indicial gust response at 
M^ - 0.9 in a superposition integral for 
arbitrary wave forms. 



Fig. 13 Indicial gust response on the 
NACA 64A0012 airfoil. 


11 



12 


NACA64A006, = 0.9 , A 1 = 1 


i-0 r 



Fig. 14 Nonlinear indioial gust response and instantaneous pressure distributions on the NACA 64A0012 

airfoil at M - 0.9 and A. = 0.017. 
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Summary and Conclusions 

The prescribed-disturbance method has proven 
to be a simple and efficient method for introduc- 
ing vortical disturbances into numerical codes for 
transonic flow analyses, including inviscid tran- 
sonic small-disturbance, Euler, and thin-layer 
Navier-Stokes formulations. The technique has 
permitted the interaction of strong concentrated 
vortices, small-amplitude distributed gusts, and 
sharp-edged gusts with stationary lifting surfaces 
to be simulated in two-dimensional transonic flow. 
The results of the different codes are in excel- 
lent qualitative agreement. However, as might be 
expected, the transonic small-disturbance calcula- 
tions can be deficient and misleading in the 
important region near the leading edge, unless 
special precautions are taken. Also, the small- 
disturbance code is less robust in coping with 
very strong vortices and with vortex paths that 
pass extremely close to the airfoil. 

Several important features of the interaction 
of a concentrated vortex with an airfoil have been 
established. Detailed examinations of results 
show that the effect of the vortex is felt pri- 
marily through the vertical velocity that it 
induces, which to first order appears as a time- 
dependent perturbation in the effective angle of 


attack, and to a lesser extent through its hori- 
zontal induced velocity. Also, unsteady lag 
effects appear to be very important, especially in 
the transonic case. Calculations show great dif- 
ferences between quasi-steady and unsteady solu- 
tions and between the results for the vortex loca- 
tions upstream and downstream of the airfoil. 
Finally, other strong nonlinear effects have been 
noted at transonic speeds. Even the qualitative 
pressure variations with respect to time and space 
differ markedly from the predictions of linear 
theory, and the vortex distorts the flow on the 
nearest surface of the body much more than the 
flow on the opposite surface. 

A vortex with an assumed, invariant structure 
has been considered thus far. In reality, close 
encounters with a body will probably alter the 
vortex core significantly and may lead to vortex 
bursting; this aspect of the problem needs to be 
examined. For most practical applications, 
including helicopter rotors, future studies need 
to be extended to include three dimensions, where 
it will be even more essential to minimize the 
level of complexity in the governing equations 
without sacrificing accuracy. 

The transonic small-disturbance formulation 
has been extended to include distributed vortical 
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disturbances superimposed on a uniform free 
stream, under the assumption that the prescribed 
disturbance remains undistorted as it convects 
past a thin airfoil. A number of unsteady air- 
foil-gust interactions have been calculated within 
this framework, and sample results illustrate the 
essential effects of periodic gusts on the tran- 
sonic flow around conventional and supercritical 
alrfoil3. Also, several examples of sharp-edged 
gusts have been studied. As in the concentrated- 
vortex interactions, the strength and unsteady 
motion of the shock wave were found to play major 
roles in the flow-field development and, conse- 
quently, in the airloads on the airfoil. 
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